Strong-Field QED and the Inverse Mellin Transform* 
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We introduce the technique of inverse Mellin transform in a problem of strong-field QED. We show 
that the moments of pair production width in a uniform background magnetic field are proportional 
to the derivatives of photon polarization function at the zero momentum. Hence, the pair-production 
width or the absorptive part of the photon polarization function is calculable from the latter by 
the inverse Mellin transform. Using the Kramers- Kronig relation, the dispersive part of photon 
polarization function can be computed as well. Therefore the analytic property of the photon 
polarization function in all energy range is obtained. We also discuss briefly the possible extensions 
of this technique to other problems. 

PACS numbers: 12.20.Ds, 11.55.Fv 



I. INTRODUCTION 

The wave function and energy quantization of a relativistic charged particle in a uniform background magnetic field 
is well understood. For a relativistic electron (positron), one solves for the Dirac equation in a background magnetic 
field. The solution is characterized by the energy levels: 

K,. m =m 2 e +p 2 z + eB(2n + l + 2s z ), (1) 

V.O ' with s z the electron spin projection along the +z direction, and the correspondent wave function 

o 

g ; *Pi, Py , Pz , Sz (r,t) = exp(-iE n t + i Py y + ip z z)F^ s ^x'), (2) 

where F^ p , fS x '^ 1S a f° ur ~ com P onen t spinor with x' = x — p y /eB. The detailed form for F* p ^ Sz can be found, for 
example, in ||]. 

In spite of our knowledge in the above energy quantization and wave function, it is often non-trivial to compute a 
physical process occurring in a background magnetic field, particularly, if there are more than one charged fermion 
involving in the reaction. A well-known example is the pair production process, 7 — » e + e , which is relevant to the 
gamma-ray attenuation in the neutron star Q , and was studied first by Toll || and Klepikov (i) independently. Both 
authors computed the pair production width by squaring the 7 — > e + e~ matrix elements using exact electron and 
positron wave functions in a background magnetic field, and summing over all available final states consistent with 
the initial photon energy. The most updated calculation using this approach was given in ||. We note that all of 
the above works assume q 2 = for the photon momentum, despite the presence of background magnetic field. To 
apply the precise photon dispersion relation in the calculation of pair-production width, one needs to study the photon 
polarization function in a background magnetic field. In fact, a consistent calculation of pair-production width require 
the knowledge of both absorptive and dispersive parts of photon polarization function. In this regard, Tsai and Erber 
H obtained the absorptive part of the one-loop photon polarization function in the asymptotic limit u> ^ 2m e and 
B B c = m 2 /e with the assumption q 2 = 0. Their result was shown Q to agree with that of Toll and Klepikov. 
However, in the above asymptotic limit, the threshold behavior of the pair-production width is completely lost. Later 
on, Shabad @ obtained the absorptive part (and the dispersive part as well) of one-loop photon polarization function 
for a general photon energy and magnetic-field strength. In that work, the threshold behavior of the pair-production 
width was worked out explicitly. We remark that Refs. |6jQ] employed Schwinger's proper-time representation for 
the electron and positron Green's functions which make up the photon polarization function. To our knowledge, 
Rcf. 0| is the first work demonstrating that the proper-time representation for the photon polarization function gives 
equivalent pair-production width to that given by squaring the 7 — > e + e~ amplitude directly. Unfortunately, the 
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manipulations of Ref. M are rather involved and the substantial details of them were given in some other unpublished 
preprints ||. It is not very clear how one can generalize the approach of Ref. |Q to other processes. 

In this report, we provide an alternative derivation of the pair-production width (or equivalently the absorptive 
part of the photon polarization function) from the proper-time representation of the photon polarization function. 
We shall also outline the procedure of obtaining the dispersive part of the polarization function, which is required for 
computing the photon index of refraction. It will be clear that our approach is very straightforward and physically 
intuitive. This report is organized as follows: In Section 2, we derive the sum rule that relates the moments of the 
pair-production width to the derivatives of the photon polarization function at the zero momentum. In Section 3, we 
apply this sum rule to the photon polarization function in an asymptotic limit, i.e., u> 3> 2m e and B <C B c . A general 
analysis valid for arbitrary lo and B is given in Section 4. We conclude in Section 5. 



II. THE SUM RULE 



We are interested in the properties of photon polarization function in the background magnetic field, II M1/ . Since a 
photon has two different polarization states in a background magnetic field, one may choose the polarization states as 
ej^ = (0, ej|) and = (0, e±) where e]| is lying on the plane spanned by the photon momentum q and the magnetic field 

vector B, while e*j_ is perpendicular to the plane. It is then convenient to define the scalar functions Ily j_ — ejj 1 
which govern the behaviors of polarization states eff ' j_ in the background magnetic field. The sum rule for II|| 
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derived by considering the following contour integral |10|: 



dco 2 Ii\\,±{^ 2 ) 
2ttT (lo 2 + w 2 )™ +1 ' 
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(3) 



where the integration contour C is shown in Fig. 1. We have implicitly assume q 2 = in the above equation in 
order to compare with previous results ||. This assumption will be relaxed later. For convenience, let us take B to 
be along +z direction, while the photon momentum q is on the XZ plane making an angle 9 to B. Hence q can be 
parameterized as q' 1 — (lo,losui9,0,locos9). We note that the integral I n may be evaluated in two different ways. 
One computes /„ either by the residue theorem or by a direct integration along the contour C with the realization 
that the contribution from the outer circle vanishes. The equivalence of two integration procedures gives rise to the 
relation: 



d" 



d(uj 2 ) n 



n. 
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. Imn |u (q; 2 ) 

(u 2 + Luir+ i 



(4) 



where Mn j_ are the threshold energies for pair productions, given by 



M,? sin 2 9 = 4m 2 , M i sin 2 9 = m 2 e 1 



; B 
'B~ c 



(5) 



Since the absorption coefficient (pair production width) is given by K\\±_ = ImII|i Eq. (W) relates the real part of 
photon polarization function to the absorption coefficient. One can further simplify Eq. (0) by taking ujq — such 
that 
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(6) 



with y\\,±_ = M 2 j_/cj 2 . One notes that the absorptive part of Hy ,j_(w 2 ) vanishes for the range < lo 2 < Af 2 ± [ pr| . 
Therefore one can effectively set the integration range of Eq. (||) as from y» j_ = to yux — oo. As mentioned before, 
the above sum rule is derived by assuming q 2 = 0. For a general q 2 , the sum rule becomes [O: 
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l,-L 



— 2n 
11,-L 



du 



U' u \\,± 



"11,-L J 



(7) 



where m 2 — 4m 2 , m\ = m 2 (l + y/\ + 2B/B C ) 2 , q 2 = oj 2 — q 2 , and u^j_ = m 2 ± /q 2 - Both of Eqs. @ and @j indicate 
that the moments of absorption coefficients are proportional to the derivatives of photon polarization functions at 
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the zero value of q 2 = to 2 — q 2 . Hence the absorption coefficients can be calculated from the derivatives of photon 
polarization functions at the zero value of q?, by the inverse Mellin transform. 




FIG. 1. The integration contour for I n and the analytic structure of IX m 
circle to infinity. 



In actual calculations, we take the radius of the 



III. THE BEHAVIORS OF IT IN THE ASYMPTOTIC LIMIT 



To study the behaviors of IIim , let us begin with the proper-time representation of photon polarization function 



ILj„ in a uniform background magnetic field p3| : 



+1 



da dvie-^itfg^-q^No 



(9\\9\\iiv ~ Q\\nQ\\u)N\\ + {q\g± l j,u - q± fl q^)N±] 
e-»'(l-« 2 )(,V- M „)}, 



(8) 



where the photon momentum has been decomposed into qjf = (a;, 0,0, q z ) and q± = (0, q x , 0). The phase cf>o and 
the functions No, N\i, and N± are given by 



t> = m, 



o 1 — V 2 2 COs(zu) — COs(z) 2 



2zsin(z) 
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with z = eBs, and 



N 

N \\ 

AT, 



cos(zv) — v cot(z) sin(zw) 
sin(z) 
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(10) 
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In the limit B «B C with q 2 = 0, we write the derivative of Him into an asymptotic series in B/B c 



d(w 2 ) n 



n,. 



2a e m 2 e fB 2 sm 2 9\" T(3n - l)r 2 (2n) 
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(11) 



where the disregarded terms are of higher order in B / ' B c 
coefficients Kii i via the inverse Mellin transform JlG]: 



From the sum rule in Eq. (1q), we obtain the absorption 
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(12) 



where a is any real number greater than 1/3; while A = (a; sin 9B/V3m e B c ) and A" = A' • (1 + 
These absorption coefficients were shown jlG] to agree with previous results by Tsai and Erber 



sin# (eB/m e ) Tim (A) with A = |(e£?/m 2 )(w/m e ) sin0, and 



l + 2B/B c )/2. 
who obtained 



7||,±(A) 



4^/3 
^rA" 



cfo(l ~v 2 y 



A' 



2/3 



A 1 



(13) 



with K 2 /3 the modified Bessel function. 

It is important to point out that both of Ku and n± are smooth functions of the photon energy. This behavior 
actually does not describe the real physical situation where new absorption peaks keep on emerging as the photon 
energy increases. One attributes this problem to the asymptotic expansion performed in Eq. (|ll|). As it is easily 
seen that, apart from the trivial mass factor m 2-2 ", the first term on the R.H.S. of Eq. (|TT| ) grows to infinity as n 
increases, no matter how small the ratio B / B c is. This then implies that the disregarded higher-order terms are in 
fact non-negligible for a sufficiently large n. As a result, the large moments of are not accurately determined by 
the first term of the above asymptotic expansion. Therefore, to obtain a correct energy dependencies of the absorption 
coefficients Ki\ ±, one should not expand in B/B c even for a small magnetic-field strength. 



IV. THE EXACT ONE-LOOP RESULT FOR U h± 

In this section, we calculate IT|| _j_ for a general photon energy and background magnetic-field strength. From 
Eq. (||), we obtain the integral representation for the scalar functions H|| _i_ , i.e., IX 1 1 (u 2 , qj_, 9) = —(a e uj 2 sin 2 0/47r)fl|| 
and U±(uj 2 ,q 2 ± ,9) = -(a e q 2 L /47r)fl i , with 



IX 1 1 = / dz dv exp[— is(f>o]N\\ , 



II_l = J dz J dv exp[-is0 o ]A ± , (14) 

where = q 2 + q 2 , and 9 is the angle between the photon momentum q and the magnetic field vector B. The 
calculations of IX 1 1 and n j_ are non-trivial. For illustrations, we will only show the details of computing one particular 
term in IX||. Let us consider the following integral 

Llf = r dz [ +1 dv exp[-^o]^$, (15) 
11 Jo J-i sin(z) 

where cos(zw)/ sin(z) is the last term in Nu. To apply the sum rule given by Eq. (0), we rotate the integration contour 
s — * —is and write ftjj 1 in an infinite series of q 2 , i.e., 
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1 (q' 2 ) n ~ l r°° 
5f (gf ,qf) = J2 K lt P A4l)(&) n -r}^ J dz z n ~ l exp[-z/?] 



(n-l)\ Jo 

x / dv (1 - v 2 ) n ~ l (exp[azv] + exp[-azv}) , (16) 



with £ = £~ „ E?=o £™=o >P =o: 9f = 9f /4m», 2?' = B/B C) = qi/4m 2 , a = p> - p + I, P = p + p' + 2m + 
1 + and 

K A ( , 2) = 2(-l) l +P+P' (2q'j)T(/ + m+l) 
imwAq_U (i-p)! J >!(f- J /)!p'!r(m + l) ' 1 ' 

Note that the various indices in the summation arise as follows: the index n comes from the photon-energy expansion, 
I arises from the binomial expansion of ((1 — v 2 )u' 2 + 2(cosh(zw) — cosh(z))q'^/zsinh(z))™, p and p' arises from 
writing (cosh(zv) — cosh(z))' as a sum of exponential functions, and m is due to expansions such as sinh - ' z — 
2 l exp[— Iz] £m=o ^m"" 1-1 ex P[ — 2mz]. After performing the integrals in Eq. ( |i~6| ) and replacing n — I + 1 by n, we 
arrive at 

n [| ^T Sm 0U W ~ °l^ K lmppi{<\±.){B) { — ) 

rfn) „ ,n n+1 1 a 2 , 
TT^rC-,— ;n + -;^). (18) 



r(n+|) v 2' 2 ' 2'/3 : 
It is easy to calculate the derivatives of Iljj 1 . Let us define -Fjj^fVi, q^, -B) = m 2n ^j d ^i^ n | g a — o ■ We have 

4 = J- f d8Ff(a,<&,B)uj', (19) 
by inverting the sum rule in Eq. (|7|). The above integral transform is evaluated using 



2 _2 , 



' / ^--IW ^(i. l±i;. + i;^, . -y'-^y. ,20, 
2tt« Jc r ( s + 2) 2 2 2 ,/l x I x z 



We arrive at 
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,Kt P A^(^-f ) + (ff) 2 ) 



^sin 2 ^ V " (21) 



a e B' 

Aq '\\" B\l(l-^) + (^) 2 

V *u ; V2 «ii 

where £" = Eto Em=o P =o- ^° understand the structure of wjj 1 , we rewrite the denominator of the above 
equation as 



with li = ((3 + a — l/B')/2 and li — [fi — a — l/B')/2. It is easy to verify that this denominator has a root at 
q 2 = m 2 (Vl + 2l\B' + \/l + 2I2B') 2 , which is precisely the threshold value for q 2 to produce an electron-positron 
pair at the Landau levels l\ and I2 respectively. This threshold behavior is also indicated by the step function in the 
numerator. 

Having worked out the q 2 dependence of /cjj 1 , we proceed to perform the summation in We first perform the 
summation over the index I in Kj^ lpp ,: 

E E tZm-Vv = £ xVr(f5 + m + ^(-^ (23) 
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where x = — 2q l2 /B' = —q 2 L /(2eB), p = m&x{p,p'}, p — min{p,p'}, and B v ^ m is the Laguerre polynomial. We then 

rewrite the left-over summation X^ p '=oSm=o as X/h=l 2^=o Em=o with t = min{ii — 1, Z2} using the relation 
li = (f3 + a — l/B')/2 = p' + m + 1 and I2 = {(3 — a — l/B')/2 = p + m. The summation over the index m can now 
be performed as follows: 

V x~ m - 1 - 1 = x _ 1 x \ (241 

^ r(r-m + i)r(f-m + i)r(m + i) r(f+i) ? v h K ] 

where I = max{Zi — 1, 12}. With the above summations over I and m, we arrive at 



D/ 00 



K^^wsin^) > v " " / (25) 

" 2 'f £S& ,/u-*S4 + «p 

witbQ 



Tf , (x) = (-l) 1+rA x rA e x ^(Aa + 1) rr A (_a.)xrx (_ a .) ( 26 ) 
V 2 v ; ^ 1 T{\ A +r A + l) Xa ' Xa ' 



where r A = I - I = | Z x - Z 2 - 1|, and = f = (Z : + ? 2 - |^i — ^2 — 1] — l)/2. 

Apart from the factor tjsin 2 6>, the absorption coefficient Kjj 1 is a function of a; = — q 2 L /(2eB) and <7y 2 where the 
latter determines the threshold behaviors of the absorption coefficient. Now that we have obtained the absorptive 
part of the polarization function , the dispersive part can be calculated using the Kramers- Kronig relation. It is 

convenient to recall the relation iljj 1 = — (a e ui 2 sin 2 0/4^)11^ (c/, 2 , q 12 )- Hence 



, r°° Imlbf(t,q' 2 ^ 



Renf(gf,q?) = - / dt 11 , (27) 

^ Ji 1 ~ 9|| 

where P stands for evaluating the principle part of the integral. We point out that the refractive index for the 
polarization state, Em 1 , can be calculated through the equation q 2 + Rell = 0. Hence 

It is important to note that the index of refraction n» also appears on the R.H.S. of the above equation, through 
relations among the quantities, q' 2 , q 12 , and the angle 9. Thus the calculation of n\\ is in principle nontrivial even 
if the polarization function ReII|| is known. Let us consider a simplified case where the magnetic-field strength is 
super-critical, i.e., B ^> B c . In this case, the quantity 1 — n| to the leading order in B c /B is given by the lowest 
Landau-level contribution to the dispersion integral in Eq. (P8|). Hence, we arrive at fll2fl 



9/ 10 n „n a e B'sin 9 , 

n|(gf ,a'l,6) - 1 = 2?rg/2 I —;===== arctan , 



<l ,\ , 2q2 



exp(-^i). (29) 



It is instructive to rewrite exp(— 2q' 2 /B') as exp(— 2lu' 2 sin #n 2 /£>') with a/ = uj/2m e . Clearly, for q 2 not very close 

to 4m 2 (the threshold energy for the pair production), one can safely set exp(— 2q' 2 /B') = 1. The index of refraction 
n 2 with this approximation then agrees with the previous result [fblf . However, as q' 2 is so close to 1 such that 

ote/ yT — gj 2 becomes greater than unity, one can no longer neglects the exponential factor cxp(— 2q' 2 /B'). In fact, 



1 The absorption coefficient is not written in a symmetrized form for saving the space. Nevertheless, the symmetrization with 
respect to h and h can be easily done as one wishes. 
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in the limit q'^ — > 1, one finds that q'j 2 — > oo, as pointed in the previous literature l^Jlq]- We note that the behaviors 
of another polarization state, ej_, can be studied in the similar manner. The details are given in Ref. p^ j. 

For consistency, it is desirable to compare our absorption coefficients with those obtained by squaring the 7 — > e + e~ 
amplitude directly^]. Our results reduce to that of Daugherty and Harding in the special limit q 2 — 0, which is an 
assumption made by these authors in their calculations. Our finding is similar to what Shabad has demonstrated in 
Ref. P] as he compared his result with that of Klepikov 0] in the same limit for q 2 = 0. Our approach differs from that 
of Refs. in that Shabad performs the calculation in the beyond-threshold energy where the algebraic manipulations 
are rather involved and precautions are required, whereas we take the advantage of inverse Mellin transform which 
permits us to calculate the polarization function near the zero longitudinal momentum with a convenient momentum 
expansion. 

V. CONCLUSION 

In conclusion, we have developed an integral-transform technique to compute the absorptive part of photon polar- 
ization function in a background magnetic field, while the dispersive part can be obtained via the Kramers- Kronig 
relation. We emphasize that this technique is rather powerful for the current problem in which there are infinite 
numbers of singularity occurring in the real positive axis of q'^ , each corresponds to an e + — e~ system in a specific 
combination of Landau levels. This technique has many other applications. To name a few, one may analyze the 
photon polarization function in a general background electromagnetic field, or study other current-current correlation 
functions under the same external condition. To our knowledge, the vector-axial vector correlation function relevant 
to weak interaction processes has not been analyzed as detailed as the fashion presented in this work. We shall report 
the results of such analysis as well as other relevant studies in a forthcoming publication |T(| . 
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